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ABSTRACT 

We present new Spitzer/IKAC observations of 55 dusty Long Period Variables (LPVs, 
48 AGB and 6 RSG stars) in the Galaxy that have different chemistry, variability 
type, and mass-loss rate. 0-rich AGB stars (including intrinsic S-type) tend to have 
redder [3.6] —[8.0] colors than carbon stars for a given [3.6] —[4.5] color due to silicate 
features increasing the flux in the 8.0 pm IRAC band. For colors including the 5.8 pm 
band, carbon stars separate into two distinct sequences, likely due to a variable pho- 
tospheric C3 feature that is only visible in relatively unobscured, low mass-loss rate 
sources. Semiregular variables tend to have smaller IR excess in [3.6] —[8.0] color than 
Miras, consistent with the hypothesis that semiregular variables lose mass discontinu- 
ously. Miras have redder colors for longer periods while semiregular variables do not. 
Galactic AGB stars follow the period-luminosity sequences found for the Magellanic 
Clouds. Mira variables fall along the fundamental pulsation sequence, while semiregu¬ 
lar variables are mostly on overtone sequences. We also derive a relationship between 
mass-loss rate and [3.6] —[8.0] color. The fits are similar in shape to those found by 
other authors for AGBs in the LMC, but discrepant in overall normalization, likely due 
to different assumptions in the models used to derive mass-loss rates. We find that IR 
colors are not unique discriminators of chemical type, suggesting caution when using 
color selection techniques to infer the chemical composition of AGB dust returned to 
the ISM. 
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1 INTRODUCTION 

The Asymptotic Giant Branch (AGB) is the last evolution¬ 
ary phase of low and intermediate mass stars (M < 8 M©), 
before the brief post-AGB phase that leads them to become 
white dwarfs. AGB stars are characterised by complex nu¬ 
cleosynthesis, high mass-loss rates, and variability. They are 
the source of a large fraction of the mass return ed by stars 
to the Interstellar Medium (ISM. ISedlmavr|[l994l) . In partic¬ 
ular, the y are believed to be the prima ry ‘dust factories’ in 
galaxies ilGehr j|l989l : lBover et al.i2012l) . even though recent 
observat ions have re-evaluated the role played by superno vae 
(see e.g. iMatsuura et al.ll201ll : lDwek fc Cherchnefilf201ll ) . It 
is clear, however, that understanding mass-loss in AGB stars 
is instrumental for modeling the chemical evolution of galax¬ 
ies. 

The composition of the material returned from AGB 
stars to the ISM depends on the chemistry of the star, which 
is determined by its initial mass and evolutionary history. 
The latter is crucially regulated by mass loss. Most stars en¬ 
ter the AGB with an intrinsic G/O abundance ratio less than 


one ( Q-rich). In stars of ~ 1.5 — 4 M© llStraniero et al.lll995l . 
119971) however, several dredge-up events following thermal 
pulses may drive the G/O ratio above unity. These stars 
become the so-called carbon stars. In stars above ~ 4 
Mq, however, carbon is destroyed by nuclear burning at 
the base of the convect i ve envelope (Hot Bottom Bur ning; 
ISmith fc LamberOl 19851 : iBoothrovd fc Sackmannlll992h pre¬ 

venting it from reaching the surface of the star where it 
could be observed. These stars never transition to the car¬ 
bon star phase. Mass limits constraining the formation of 
carbon stars depend on metallicity; low metallicity envi¬ 
ronments (e.g. the Magellanic Glouds) form carbon stars 
at lower masses. Intrinsic S-type stars share many charac¬ 
teristics with carbon stars such as similar masses and a rich 
s-element chemistry from many dredge-up events, but they 
lack sufficient carbon to drive the C/O ratio above one. 

Stars with masses ~ 8 — 25 Mq do not enter the 
AGB phase, but become Red Supergiants (RSGs) on the 
way to ending their lives as Type II s upernovae (e.g. 
iMassev fc Olsen! I2OO3I : Ivan Dvk et al.l FooJ ) . RSGs do not 
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Figure 1. Images of O-rich AGB star CZ Ser in each of the IRAC bands before PSF subtraction (top panels) and after PSF subtraction 
(bottom panels). The horizontal and vertical bands are un-subtracted electronic artefacts that are ignored in the fit. The PSF was fit to 
unsaturated diffraction spikes and rings, outside the saturation radius (~ 20 " for CZ Ser). The field of view of each panel is ~ 5^ 


experience third dredge up episodes and a s such they are 
chara cterised by an O-rich chemistry (e.g. IVerhoelst et al.l 
l2009t) . More luminous than AGB stars, but with similar 
spectral types, they are also subject to variability and mass 
loss. Their dust production rates are comparable to those of 
AGB stars but their overall yield to the ISM is lower, due to 
RSGs being mor e rare than their lower mass counterparts 
in the AGB le.g. iBover et al.ll2012l l. 

AGB and RSG stars lie close to an instability region of 
the Hertzsprung-Russell diagram and as such they exhibit 
long period variability of Mira, semiregular, and irregular 
type with periods on the order of hundreds to thousands 
of days. Period-luminosity diagrams of long period variables 
(LPVs) in the Magellanic clouds (LMG and SMG) reveal 
sequences that appear to correspond to the variable pulsa¬ 
tion mode, bina rity and other as yet unknown characteris¬ 
tics dWoodlbOld ). These diagrams can be a powerful tool for 
studying the pulsation mode for a large population of vari¬ 
able stars. This is especially true when the available light 
curves are sparse and poorly sampled, making it difficult 
to detect additional periodicity (e.g. overtone pulsations) 
in the amplitude modulations of the light curve. The con¬ 
struction of similar diagrams for Galactic AGB stars would 
be extremely useful, as they would illuminate the role of 
metallicity in determining the pulsation properties of these 
variables. Efforts in this sense depend critically on accu- 
rate and precise distan ce measurements for Galactic LPVs 
dWhitelock et al]l2008l l. 


Pulsations and ma ss loss are inex tricably linked in AGB 
and RSG stars (see e.g. I Willsonll2000ll since radial pulsations 
aide in the formation of the dust-driven winds responsible for 
the characteristic high mass-loss r ates of AGB stars (10~^ 


Ivan Loon et al.lll99£ 

, as compared with 10 ® Mjri yr ^ to 

10“'^ Mq yr"\ e.g. 

Mauas et al.l I 2 OO 6 I: iDuoree et al.l l2009l 


for Red Giant Branch stars). The dusty wind that results 


from these mass-loss processes, before elements synthesized 
in AGB stars are released into the ISM, leads to the for¬ 
mation of dusty cocoons enshrouding the star. These cir- 
cumstellar envelopes obscure optical radiation, making ob¬ 
servation of mass-losing AGB stars difficult at visible wave¬ 
lengths. However, thermal radiation from the dust makes 
AGB stars extremely luminous in the infrared (IR ). The In- 
fraRed Array Gamera (IRAC, iFazio et al.l 200j) on board 
the Spitzer Space Telescope dWerner et al. 120041 ') operates 
in the mid-IR and is an ideal instrument for characterizing 
AGB stars according to their chemical composition, mass- 
loss rate, and variability class. IRAC’s four channel imag¬ 
ing during the Spitzer cryogenic mission was especially well- 
suited to observe the prominent dust and molecular spectral 
features in AGB stars - for example, H 2 O, SiO, CO 2 , and 
CO among others for O-rich stars and C 2 H 2 , HCN, CS, and 
C 3 for carbon stars. 


Due to their intrinsic brightness and red colors, AGB 
stars are easily detected by IRAC and have been found in 
numerous IRAC surve ys of the Milky Way (e.g . GLIMPSE, 
iBeniamin et al.l I2OO3I : iGhurchwell et aO 
Group galaxies (e.g. SA GE. lMeixner et al 


20031 and Local 


20061: iBlum et al.| 


l2006l : lBover et ^I^Olll l. However, there have been no stud¬ 
ies specifically designed to characterise in the IRAG bands 
a sample of Galactic AGB stars with well determined chem¬ 
ical type, variability, and mass-loss rate. To fill this gap, we 
designed a program to observe a nearby sample of Galactic 
AGB stars with IRAG in order to better identify AGB stars 
in IRAG surveys, to facilitate Galactic population synthesis 
studies, and to better understand the chemical evolution of 
the diffuse matter in the ISM. 


We have observed 48 Galactic AGB stars, representing 
each of the main types of AGB stars - O-rich, S-type, car¬ 
bon stars, Mira, semiregular, and irregular variables as well 
as 6 RSGs of semiregular and irregular variability type. We 
outline our target selection in Section 2. Section 3 describes 
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the data reduction pipeline and how the photometry was de¬ 
rived. We present the IRAC colors of AGB stars in Section 
4 and discuss period, magnitude and color relations in Sec¬ 
tion 5. Relationships with the mass-loss rate are explored in 
Section 6. We summarise our conclusions in Section 7. 


the need for single epoch photometry with better precision 
than 2MASS. 

Given the variability of the sources, we requested two 
epochs (six months apart, as constrained by the Spitzer vis¬ 
ibility windows) for each target, in order to check for vari¬ 
ations in the IRAG photometry on time-scales of several 
months. 


2 TARGET SELECTION 


Our target list was selected from a number of Galactic 
AGB star catalogs available in t he literature (|LoirB_et_aJj 
I 993 I: iKerschbaum fc Hronlll996l: lAdelman fc Maher! 19981 : 


Heras fc Honv 20051 : iGuandalini et al.l 20061 . and references 
therein) with the intention of representing all main types of 
AGB stars. Two main constraints limited our choice of tar¬ 
gets: (1) the availability of distance estimates and (2) a K 
band magnitude > 0, to prevent excessive saturation with 
IRAG. While our target list is not a statistically signihcant 
sample of the Galactic AGB population, the sample is suf- 
hciently large and diverse to allow a general study of the 
spectral energy distribution (SED) of Galactic AGB stars 
in the Spitzer/IKAC bands, and its overall dependence on 
their chemical, physical, and variability characteristics. 

Distances for our target stars were derived either 
from interfer o metri c observations (e.g. Ivan Belle et al ]l2n02l: 
IZhang et al.1 l2012[ l , m odels of radio emission and/or 
bolometric luminosity ( Loupe£_alJ_ W9^;_ 01ofesone£_^ 
2 OO 2 I: IGuandalini et al.l 20061 : IGuandalini fc Bussd 1200^ 


Schdier et al.H2013l) . a nd, in the absence of either of those, 


astrometric methods dvan LeeuwenI [ioo^ . adopting the lat¬ 
est corrections of Hipparcos distances). The distance of our 
selected AGB targets varies from 0.14 to 1.85 kpc (with one 
single carbon star that has a distance of 4.21 kpc). The dis¬ 
tances to previously observed supergiants have been derived 
with different methods found in the literature (see notes on 
Table [1]) . 

The complete sample is listed in Table [U and is di¬ 
vided in three main categories: 22 O-rich AGB stars (M III 
spectral type), 7 intrinsic S stars, 19 carbon stars, and 6 
supergiants. In each category we have a similar number 
of Mira, semiregular and irregular variables, with periods 
ranging from 50 to 822 days. Information about the pe¬ 
riod and variability type of the targets was o btained from 
the G eneral Gatalog of Variable Stars (GGVS. ISamus et al.l 
l2012h . Estimated mass-loss rates were obtained from radio 
observations of CO or HCN in the outer envelope, or fit¬ 
ting detailed radiative transfer models of the circumstel- 
lar emission to the infrared SED (see ILoup et al.1 1 19931 : 
lOlofsson et al.l[2002l : IGuandalini et al]|2006l . and references 
therein). The mass-loss rates of our target stars range from 
10“® to 10~^ Mq yr“^ in each category. Mass-loss rates 
for the supergiants were estimated by various methods (see 
notes in Table [!}. Near-IR photometry w as obtained (for 
rand om epochs) from the 2MASS catalog dSkrutskie et al.l 
I 2 OO 6 I I. While higher quality near-IR photometry exists for 
some of the stars in our sample, we chose to adopt 2MASS 
photometry for uniformity across the sample. The uncer¬ 
tainty of the 2MASS magnitudes are already smaller than 
the amplitude of the infrared variability of our LPVs. In 
absence of complete light-curves for all stars in the sample 
from which to calculate average magnitudes, this obviates 


3 OBSERVATIONS AND DATA REDUCTION 

We restricted our target list to nearby AGB stars with 
estimates of their distance. Given the intrinsic luminosity 
(~ 10^ Le ) and red colors of AGB stars, these targets sat¬ 
urate the IRAC detectors even at the shortest IRAC subar¬ 
ray frame times. While saturated images are not suitable for 
aperture photometry, reliable photometry can be recovered 
with Point Spread Function (PSF) fitting, as long as at least 
part of the PSF is not saturated. To compromise between 
the need to limit the amount of saturation in our images and, 
at the same time, fill enough of the IRAC field-of-view with 
high S/N unsaturated portions of the PSF, we adopted the 
2 sec full-frame IRAC Astronomical Observation Template. 
To allow for efficient removal of outliers (bad pixels and cos¬ 
mic rays), and sufficient spatial sampling of the unsaturated 
PSF optical features (diffraction spikes and rings), we ob¬ 
served each target using the 5-point small scale Gaussian 
dither pattern. 

The data were acquired between 2006 June 2 and 2008 
June 20 as part of the IRAC Guaranteed Time, with PID 
30411. Starting with the Basic Calibrated Data (BCD) pro¬ 
duced by the Spitzer Science Center (SSC) pipeline version 
S14.4.0, we generated a mosaiced image for e ach source, us¬ 
ing o ur own post-BCD software IRACproc dSchuster et al.l 
I 2 OO 6 II . IRACproc is an add-on to SSC’s MOPEX mosaicing 
software, that applies more sophisticated outlier rejection 
criteria. With the optimised parameters from IRACproc, 
MOPEX defines a fiducial image frame that contains the 
coordinates for the constituent frames, removes transients 
using temporal outlier rejection and then interpolates the 
frames to produce a mosaic image. The final mosaic is scaled 
to 0.863”/pixel which is l/%/2 of the IRAC pixel scale (half 
of the IRAC native pixel area), to provide ideal (Nyquist) 
sampling even at the shortest IRAC wavelengths. 

Table [3 lists the measured IRAC Vega magnitudes for 
all sources. The adopted PSF fitting technique, developed 
sp ecifically for he a vily s aturated IRAC images, is described 
in iMarengo et all 1 I 2 OO 9 I '). As an example of this procedure, 
Figure [ 1 ] shows the O-rich AGB star CZ Ser, shown before 
and after the PSF is fit and subtracted from the target. We 
used a high dynamic range image of the IRAC PSfQ in each 
band, made from a combination of individual images of a set 
of stars with different brightness (Sirius, Vega, Fomalhaut, 
e Eridani, e Indi, and the IRAC calibrator BD+68 1022). 
The PSF intensity is scaled to match the actual IRAC ob¬ 
servations of Vega in each band, providing an absolute pho¬ 
tometric reference, and is super-sampled over a grid with 
0.24"/pixel. By htting the unsaturated parts (diffraction 


^ available at the IRSA website: 

http: //irsa. ipac.caltech.edu/data/SPITZER/docs/irac/calibrationfiles/psfprf/ 
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[4.5] - [8.0] [5.8] - [8.0] 

Figure 2. Color-Color diagrams (average of the two epochs) of the observed AGB stars. O-rich stars are marked with d i amon ds, carbon 
stars with triangles, intrinsic S-type stars with squares and supergiants with crosses. Model tracks from ICroenewe^ 3 il2onfit) for three 
spectral types are overplotted. Carbon star models are plotted with dotted lines, and models with an MO III and MIO III central star 
and silicate dust are plotted with solid and dashed lines, respectively. 


spikes and PSF ring and “tails”) of each target star image, 
we have measured the flux ratio between the stars in our 
sample and Vega. The typical uncertainty of this procedure 
is 2—5% (constrained by comparing the amplitude of the 
PSF subtraction residuals with the background and source 
photon noise), and is comparable to the typical absolute 
aperture photometry of unsaturated stars with IRAC. The 
photometry of the six supe rgiants from the literature ob¬ 
tained from ISchusterl |2003) use the same PSF-fitting tech¬ 
nique. 


4 IRAC COLORS OF GALACTIC AGB STARS 

Ground-based mid-IR colors have been used extensively to 
characterise the physical and chemical prope rties of the cir- 
cums tellar envelopes of evolved stars (see e.g. lMarengo et HI 
I 1999 I and references therein). The IRAC bands offer a sim¬ 
ilar opportunity, as they span a wavelength range rich in 
strong molecular and dust features, and they are free from 
variable telluric absorption lines. These features correlate 
with the chemical signature of the circumstellar environ¬ 
ment, dust mineralogy, and abundance in the stellar wind. 
Note that in the IR, interstellar extinction is dramatically 


reduced with respect to the visible; as a consequence the 
mid-IR and near-IR colors of our sources are mainly deter¬ 
mined by photospheric and circumstellar features. 

Figure [ 2 ] shows a number of IRAC color-color diagrams 
for all target sources. Each source is plotted once, using 
the average color of the two observed epochs. The sources 
tend to be organized on a sequence of increasing excess 
col ors, matching t he ov erall distribution of colors found 
by iMarengo et al.l (l2007fl using synthetic IRAC photome¬ 
try derived from Infrared Spac e Observatory Short Wave- 
length Sp ectrometer (ISO S WS, [Valentiin et al.ll 1993) spec¬ 
tra (from isloan et aDl2003ll . The largest dispersion is in the 
[3.61 —[8.01 color (panel a). Radiative transfer models (see 
e.g. iGroenewegenI I^OOal show that this sequence is related 
to increasing amounts of circumstellar dust, indicating larger 
mass-loss rates (assuming wind velocities, dust opacities, 
and a uniform dust-to-gas mass ratio across the sample). 
We explore quantitatively the correlation between mass-loss 
rates and IRAG colors in section [6| 

Stars extending along the infrared excess sequence tend 
to separate according to their circumstellar chemistry. O- 
rich and intrinsic S-type AGB stars, as well as the super¬ 
giants, tend to have a redder [3.6] —[8.0] color (for the same 
[3.6] —[4.5] color), due to the presence of silicate features dra- 
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[5.8] - [8.0] 

Figure 3. Color-color diagram of our carbon stars with different 
symbols for increasing mass-loss rate. For each star we show the 
color for both epochs, connected by a line. Stars with the highest 
mass-loss rates are on the upper branch, with redder [4.5] —[5.8] 
color than the lower branch. Stars with lower mass-loss rates are 
on the lower track, indicating depressed flux at 5.8 pm. 

matically increasing their flux in the 8.0 pm band. This sep¬ 
aration is only effective for sources with moderate infrared 
excess, e. g. [3.6] — [8.01 >J_mag. 

The lOroenewegenl (l200fil 'l models we selected miss the 
location of the most extreme O-rich AGB stars in all of the 
color-color diagrams. This may be due to the choice of dust 
chemistry in the model. In Figure O we plot models us¬ 
ing a mixed chemistry - 60% aluminu m oxide (amorphous 
porous AI 2 O 3 ) and 40% silicates (see iGroenewegenI bood . 
for details). More recent papers have shown that the choice 
of optical constants aff ects the predicted flux of the silicate 
feature near 10 pm fsee lGroenewegen et al.ir2009l l. This sug¬ 
gests that the particular silicate opacities employed in the 
model may lead to the inaccurate determination of the dust 
mass-loss rate. 

Plots including colors with the 5.8 pm band (panels b 
and d) instead show an unexpected separation of carbon 
stars into two distinct branches. Both branches show simi¬ 
lar [5.8] —[8.0] excess, but the upper branch continues to get 
progressively redder in both ]4.5] —[5.8] and [3.6] —[4.5] col¬ 
ors, whereas stars in the lower branch tend to get bluer in 
the [4.5] —[5.8] color. This divergence suggests the presence 
of a broad absorption feature at 5.8 pm. Looking at just 
the carbon stars in [5.8] —[8.0] versus [4.5] —[5.8] color and 
using different symbols to indicate increasing mass-loss rate 
(Figure|3}, we see that the separation is related to the mass- 
loss rate listed in Table [T] Stars with higher mass-loss rates 
are on the upper, red branch, while stars with lower mass- 
loss rates populate the lower, blue track. This dichotomy 
suggests that the feature responsible for the 5.8 pm absorp¬ 
tion is likely photospheric, being masked by thicker dusty 
envelopes in sources with higher mass-loss rates. 

To identify this feature, and to match the observed 
trends in the IRAG colors wi th specific spectra l features, 
we retrieved archival spectra dSloan et al.l I 2 OO 3 I 1 for all of 
our sources that were observed with ISO SWS. We have 





Figure 4. ISO SWS spectra of representative stars in our sample. 
Top: — carbon stars; Center: - O-rich AGB stars; Bottom: - two 
intrinsic S-stars and a red supergiant (AA Cyg). The IRAC pass- 
bands are plotted below the spectra and the main p hotospheric 
and circumstellar features are ind i cated (following, e.g.lAoki et al.l 
I 1999 I : iNoriega-Crespo et al.ll2004l : IZhang fc Jiandi200sll . 


also compared our colors with the synthetic photometry of 
12 2 AGB and RSG stars of different chemical type, derived 
bv iMarengo et all ll2007l] by convolving available ISO/SWS 
spectra with the IRAC passbands. Some of these sources 
have been observed by ISO at multiple epochs allowing a 
direct test of the effect of variable spectral features on the 
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IRAC colors. Representative examples are shown in Fig¬ 
ure [I] The spectrum of carbon stars along the blue branch 
(e.g. VX And in the top panel) exhibits a strong C 3 ab¬ 
sorption feature centred at 5 pm. This feature is not present 
in the spectra of star s in the redder bran ch (e.g. V CrB 
and T Ind). Based on ISloan fc Pric3 (Il998f ). C 3 is assumed 
to be developing in the atmospheres of dust poor carbon 
stars. This feature disappears in stars with thicker circum- 
stellar envelopes, either because it is Hlled by the continuum 
dust emission, or becau se the molecule makin g the feature 
is depleted. As noted in iMarengo et all |20o3), this feature 
appears to be transient, as some sources observed with ISO 
SWS in multiple epochs do not show it in all spectra. This 
suggests that the feature variability is related to changes 
in the C 3 abundances in the stellar atmosphere as the star 
pulsates, rather than changes in the dust content of the cir- 
cumstellar envelope (which is unlikely to show large scale 
variations on the short time-scales of the ISO repeated obser¬ 
vations, altho ugh such variations are found on longer time- 
scales, see e.g. IWhitelock et al.ll200al . 

The spectra of carbon stars present several other molec¬ 
ular features in the IRAC passbands (CO, C 2 H 2 , HCN 
and CS). C 2 H 2 falls mostly outside the 3.6 pm band, and 
does not contribute significantly to the [3.6] —[4.5] color of 
the carbon stars. CO absorption in the 4.5 pm band more 
likely contributes to the negative (as low as ~ — 0.2 mag) 
[3.6] —[4.5] color of the carbon stars with low overall infrared 
excess. All these molecular absorption features are filled by 
dust continuum emission for the redder sources, explaining 
the general trend of increasing [3.6] —[8.0] color for stars with 
larger circumstellar dust content. 

The mid-IR spectra of the O-rich sources (middle 
panel), intrinsic S-stars and a red supergiant (bottom panel) 
are also rich in several molecular features (including CO, 
H 2 O, SiO, CO 2 and SO 2 ). Of these, the features that have 
the largest effect on the IRAC colors are SiO, CO 2 , and CO 
bands that can severely depress the 4.5 pm flux leading to 
blue [3.6] —[4.5] colors (as low as ~ —0.4 mag) for the sources 
with less overall infrared excess. Sources with high dust con¬ 
tent show the prominent 10 pm silicate feature in emission. 
As mentioned before, this feature is partially captured by 
the 8.0 pm band, causing the [3.6] —[8.0] color to be redder 
for higher mass-loss rates. There are no significant differ¬ 
ences between the colors of the individual types of O-rich 
sources, as evidenced by the fact that M-type AGB stars, 
intrinsic S-stars and red supergiants all trace the same color 
sequence. 

A number of classification schemes for the chemistry of 
AGB stars rely on the co mbination of IRAG with near-IR 
colors. Figures 7 and 8 in iBover et al.l (1201 ill , for example, 
show how O- and C-rich sources in the SMC and LMC seg¬ 
regate according to their J—[8.0] color. With this scheme, 
O-rich AGB stars have J — [8.0] < 1.4 mag and G-rich AGB 
stars have J — [8.0] > 1.9 mag, with the so-called “anoma¬ 
lous” O-rich population having intermediate color. We do 
not see this segregation for Galactic AGB stars, as shown 
in our Figured where G-rich and O-rich AGB stars overlap 
above J — [8.0] > 2 mag and intrinsic S-type stars (proposed 
as one possible explanation for the LMG/SMG anomalous 
O-rich AGB class) overlap with the range of the anomalous 
O-rich AGB stars (see box in Figure [5]). 

The fact that we have several O-rich AGB stars that 
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Figure 5. Color-magnitude diagram using near-IR photometry 
from 2MASS and Spitzer 8 pm. The tr apezial box f r ames anoma- 
lous O-rich AGB stars as identified bv iBover et al.l ll201lh . 


have significantly redder colors than most O-rich AGB stars 
in the Magellanic Clouds may be the consequence of the 
different evolutionary paths of the two AGB populations in 
galaxies with different metallicity. AGB stars are observed 
to have redder colours at higher metallic ities because the gas 
to dust ratio increases (Ivan Loonll 200 (]li . We note that the 
few O-rich AGB stars in the Magellanic Clouds that have 
redder colours tha n the Galactic AGB stars in our sampl e 
are OH/IR stars (IWood et al.lll99^ : iMarshall et al. 2004h 
As noticed by many au thors (see, e.g. Ferrarotti fc Gail 
I 2 OO 6 I : IVentura et al.ll 2012 li a low metallicity environment fa¬ 
vors the creation of carbon stars early in the evolution of 
an AGB star, when the mass-loss rate (and hence infrared 
excess) tends to be smaller. According to this hypothesis, 
LMG/SMG O-rich AGB stars are restricted to low J—[8.0] 
excess, with only C-rich stars having redder J—[8.0] colors 
due to their greater mass loss. This restriction would not 
be present in Galactic O-rich AGB stars, requiring a higher 
number of thermal pulses to transition to C-rich chemistry. 
Alternatively (though not in our sample), the neat color sep¬ 
aration observed in the LMG/SMG C- and O-rich AGB stars 
could be explained by the selection method used to identify 
the chemical type of the source. LMG/SMG AGB catalogs 
are typically base d on near-IR colors (see e.g. Figure 5 in 
iBover et al. [20]^ while the chemistry of our targets has 
been identified spectroscopically. This second possibility is 
supported by the large range in infrared colors exhibited by 
spectros copically identified C- and O-rich LMC AGB stars 
(see e.g. iMatsuura et al.ll2009l L Finally, this may be a se¬ 
lection bias in our sample as we have selected sources with 
measurable mass-loss rates, and hence a larger infrared ex¬ 
cess. 


The two reddest O-rich AGB stars in our sample are 
KU And and V1300 Aql. Both have red ]3.6] —[8.0] color 
and VI300 Aql is the reddest in the near-IR (consistent with 
its fainter J-band magnitude indicative of higher circumstel¬ 
lar extinction, see Table [1]). Several (5) of our O-rich AGB 
stars have [5.8] —]8.0] vs. [3.6] —]4.5] color sufficiently red that 
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Figure 6. Color-magnitude diagram using near-IR photometry 
from 2MASS. Line s correspond to the near-IR color cuts used by 
iBover et alJ |20^). The red striped box denotes the RGB. 

they w ould be considered extreme AGB stars bv iBover et al.l 
(l201ll) . In fact, all of our 0-rich AGB stars are redder than 
the J — K color cut iBover et al.l |20l3) use to isolate O-rich 
AGB stars in the LMC (see Figure [ 6 |. 

To test whether the redder colors of our O-rich AGB 
stars are representative of the Galactic AGB population, or 
rather are the consequence of selection effects, we extracted 
the 2MASS p hotometry of all Mi ra and semiregular variables 
in the GCVS JSamus et ahlboi^ l with known chemical type. 
The median J — K color (0.85 ma g) of this larger sam ple of 
O-rich AGB stars falls within the IBover et al.l (I2OIII) color 
cut, suggesting that the red colors of the O-rich AGB stars 
reflect our choice of O-rich AGB stars with measured mass- 
loss rates, and thus greater IR excess (see Figure[7|. It should 
be noted, however, that a large fraction of the GCVS carbon 
sta rs have a. J — K c olor falling outsid e the boundaries o f 
the IBover et al.l | 20 l 3 ) G-rich color cut. ljavadi et al] (I2OIII ) 
also find that carbon stars are less red at similar metallicities 
in M33. We find that 78% of the G-rich AGB stars in the 
Galaxy are bl uer than J — K ^ 1.2 mag (the color selection 
criteria from IBover et al ] l 201 ll. for K~ —8 mag, see Fig¬ 
ure El). Based on this, the majority of C-rich sources would 
likely be misclassified as O-rich sources. However, sources 
in the GGVS tend to be optically selected, and as such are 
biased toward low circumstellar extinction (small mass-loss 
rates), and may sample less-evolved, non-dusty AGB stars 
that behave v ery differently than their more evolved, dustier 
counterparts JBover et aklboilh . 

4.1 Source Variability 

The availability of two separate epochs, while not sufficient 
to reconstruct a full light-curve, is enough to estimate the av¬ 
erage variations of LPVs in the IRAC bands. Figure [ 8 ] shows 
the histograms of the average magnitude change between 
epochs for all sources in our sample. The RMS variations 
are 0.29, 0.28, 0.27, and 0.26 mag at 3.6, 4.5, 5.8 and 8.0 pm 
respectively. This is much larger than the RMS variations 



-10 12 3 

J-K 


Figure 7. Histogram of the J — K color of all AGB stars of 
known chemical and variability type from the GGVS. J and K 
photometry from 2MASS. The shaded region indicates the J — K 
color range of the O-rich AGB stars in the IBover et al.l ll201lh 
color selection. 


expected for a sample of non-variable stars (< 0.008 mag for 
all bands), based on the average photometric uncertainty. As 
expected, the IRAC bands are affected by the variability of 
the sources to a much smaller extent than the optical bands 
(whe re Mira light-curve s can have amplitudes as large as 11 
mag, ISamus et al.l 12012 ). with smaller variations for longer 
wavelengths. The smallest changes are observed for sources 
that are classified as semiregular or irregular variables (0.14, 
0 . 12 , 0 . 11 , and 0.10 mag in each band), as expected because 
of their shorter periods and smaller pulsation amplitudes 
with respect to Miras (0.46, 0.45, 0.42, and 0.41 mag re¬ 
spectively). 

The IRAC color changes between epochs are generally 
small (RMS ~ 0.08 mag), resulting in a small scatter in the 
color sequences shown in Figure [2] The largest variation is 
found in the [5.8] —[8.0] color for the carbon stars in the blue 
branch (RMS ~ 0.11 mag), affected by the variability of the 
C 3 feature (see figure [3]). One carbon star (RS Cyg) also 
shows a first epoch [3.6] —[8.0] color significantly bluer (by 
0.31 mag) than all other carbon stars. Unfortunately, ISO 
SWS spectra are not available for this source. However, we 
found that our second epoch IRAC photometry closely fits 
the spectrum of another carbon star, T Ind. The first epoch 
IRAC photometry of RS Cyg is characterised by a much 
higher 3.6 pm flux than the second epoch, or the T Ind 
spectrum. We suspect that this fluctuation is due to unusu¬ 
ally strong, variable C 2 H 2 absorption, common in many of 
the ISO SWS spectra of the other carbon stars in our sam¬ 
ple. L and M band spectra of RS Cyg taken at multiple 
phases would be required to confirm this hypothesis, and 
the existence of such broad variations in the strength of this 
absorption feature. 
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Figure 8. Histogram of the magnitude variations between the two epochs of observation. White boxes show variations of the full AGB 
sample, with Mira variables indicated with gray boxes, semiregular and irregulars indicated with striped boxes. 


5 PERIOD, MAGNITUDE AND COLOR 

RELATIONS 

LMC and SMC LPVs tend to organize themselves in a 
series of parallel sequences i n optical and infrared period- 
luminosity diagrams (see e.g. IWoodlli^Old c lRiebel et al.ll201Cll 
and references therein). These sequences depend on the pul¬ 
sation mode, binarity, and other not yet identified parame¬ 
ters (giving rise to poorly understood characteristics, like a 
long secondary period observed in some stars). The so-called 
C sequence is populated by fundamental mode pulsators. 
The brightest of the C sequence variables tend to be Miras, 
while the lower end of the sequence (straddling the RGB 
tip) is mainly inhabited by semiregulars. Sequence C’ (to 
the left of sequence C by A[logP] ~ —0.1), is instead pop¬ 
ulated by first overtone semiregular variables. LPVs with 
higher overtone modes are organized in separate sequences 
(A and B) with even shorter periods. Stars with the myste¬ 
rious secondary long period are found in sequence D, shifted 
by A[logP] ~ -1-0.3 from the C sequence. 

Figure [5] shows the period-luminosity diagram calcu¬ 
lated for the 3.6 pm absolute magnitude (averaged between 
the two epochs) for our sample of AGB stars. Diagrams for 
the other t hree bands are sim ilar. The sources classified in 
the GCVS dSamus et al. Iffl) as Mira and semiregulars are 
indicated by different symbols. The size of the symbols is 
proportional to the [3.6] —[8.0] color (as a proxy for their 
mass-loss rate). The shaded bands for the individual se¬ 
quences are derived from the approximate distribution of 


LMC LPVs in iRiebel et al.l (l2010l l where the [3.6] magni¬ 
tude is taken as a proxy for the luminosity, and converted 
to ab solute magnitudes using the LMC d istance modulus of 


to ab solute magnitudes using tne Ll 
18.5 ll Freedman et al.ll200ll : I^Mlmi] 


are similar to those found bv lWoodI ll201Cl) using the K band 


20n). These sequences 


magnitude as a proxy for luminosity, however they have a 
slightly shallower slope (larger period for a given luminos¬ 
ity)- 


The location of our sources on the diagram is listed 
in Table [3] and is largely as expected, based on AGB 
stars in the Magel lanic Clouds, confirming results found by 
Glass et al.l (l2009ll us ing synthetic IRAC photometry from 
Marengo et al. IpOOTfl . Miras (squares) are all within or be¬ 
low the sequence C (fundamental pulsators), with the excep¬ 
tion of the carbon star R For, which is in the first overtone 
sequence C’. A few of the Miras appear to be significantly 
under-luminous in the 3.6 pm band, up to 3 mag below the 
C sequence for their period, placing them in or near the 
secondary long period sequence D. Most of these stars (e.g. 
KU And and V1300 Aql) have large infrared excess due to 
their high mass-loss rate; their low 3.6 pm brightness may 
be caused by the large extinction due to their thick circum- 
stellar envelope. Alternately, they may have increased their 
period as the stellar structure responds to the mass lost. The 
semiregular variables, as expected, are mostly located in the 
overtone sequences or the lower part of sequence C. Only 
four semiregulars (T Ind, VX And, BN Mon and RS Cyg, 
all carbon stars) are found in the area of the C sequence oc- 
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Figure 9. 3.6 |rm absolute magnitude (average of the two 
epochs) versus the logarithm of the period in days for all Miras 
(squares), semiregular AGB variables (triangles), and supergiants 
(diamonds) with a period available in the literature. Shaded and 
str iped boxes in d icate t he approximate distribution of LMC LPVs 
in iRiebel et al.l ll2ni0h . Symbol sizes correspond to [3.6] —[8.0] 
color (taken as a proxy for their mass-loss rate). 

cupied by fundamental mode Miras. ET Vir is the faintest 
source in the plot, lying at the very bottom of sequence C. 
At a Hipparcos-determined distance of 0.14 kpc, ET Vir is 
also the nearest source in the sample. Given the large un¬ 
certainty on the parallax (> 30%), we might expect that 
the luminosity has been underestimated due to an unreli¬ 
able distance estimate. However, the Hipparcos value agrees 
well with the distance modulus calcu lated using th e width of 
chromospheric Ca ll emission lines llWilsorJll976l ) to deter¬ 
mine the absolute magnitude, making it unclear why ET Vir 
is more than a magnitude fainter than the other sources in 
the sample. One possibility is that we observed ET Vir with 
lower luminosity following a thermal pulse. All semiregular 
AGB variables in our sample tend to have small infrared 
excess, so it is unlikely that their position in the diagram 
is significantly affected by extinction at 3.6 pm. The only 
semiregular variables with large excess are four supergiants 
(XX Per, W Per, NML Cyg, and VX Sgr), located at the 
top of sequence C’ as expected. 

The general distribution of our Galactic LPVs in the 
period-luminosity diagram is characterised by a significantly 
larger spread than the variables in the Magellanic Clouds . 
As found by previous authors (see e.g. iTabur et al.l 120101) . 
this is primarily due to the larger uncertainty in the dis¬ 
tance of Galactic LPVs, and possibly in the larger spread 
in their metal abundance, which also affects their period- 
luminosity relation. Despite these difficulties, the period- 
luminosity plot using the IRAC 3.6 pm band (where ex¬ 
tinction is minimal, and infrared excess tends to be lower 
than at longer wavelengths) as a proxy for luminosity re¬ 
mains an important tool to study the pulsation mode of 
these variables, most of which only possess sparse and in¬ 
complete light-curves. 

While period-magnitude relations are diagnostics of the 
pulsation mechanisms of several classes of variables, the in- 


Figure 10. IRAC [3.6] —[8.0] color (average of the two epochs) 
versus the logarithm of the period in days. A fit to the semiregular 
variables (triangles) shows a strong correlation with period (dot¬ 
ted line). A fit to the Miras (squares) is less well-defined (dashed 
line). 

frared colors of AGB stars are dominated by dust emission 
and molecular absorption. The relationship between infrared 
colors and period, however, provides an excellent diagnostic 
to investigate the dependence of the circumstellar chemistry 
and mass-loss rate on the pulsation period. The [3.6] —[8.0] 
color, which is strongly affected by dust emission, is a good 
candidate for this analysis. 

The results are shown in figure[T0l Semiregular variables 
are distributed along a tight (RMS ~ 0.12 mag) period-color 
relation, statistically consistent with a zero slope, confined 
within a narrow color range (A([3.6] —[8.0]) < 0.5 mag). 
Note that all semiregulars, including the ones with larger 
period that appear to be on the fundamental mode sequence 
C, are part of the same period-color relation. The Mira vari¬ 
ables are instead spread over a much larger range of colors 
(A([3.6] —[8.0]) Rs 3 mag), with longer period sources (on 
the D sequence) having the largest infrared excess. The two 
linear best-fitting relations are: 

Semiregular 

[3.6] - [8.0] = (-0.11 ± 0.10) -b (0.17 ± 0.05) log(P) (1) 
RMS = 0.12 mag 

Mira 

]3.6] - [8.0] = (-7.93 ± 0.41) + (3.45 ± 0.16) log(P) (2) 
RMS = 0.42 mag 

where the period is in days. 

One possible explanation for this dichotomy is that the 
two groups of variables are characterised by different dust 
composition, affecting their excess in the IRAG bands. One 
problem with this interpretation is that there seems to be no 
difference between 0-rich (M and S-type) and G-rich AGB 
stars, despite the very different source of the infrared excess 
(mainly continuum from amorphous carbon dust for carbon 
stars, and the 10 pm silicate feature in emission for M and 
S-type AGB stars). Another possibility is a difference in the 
amount of circumstellar dust responsible for the IR excess. 
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i.e. a different mass-loss rate. While the median mass-loss 
rate of the Miras in our sample is indeed higher than for the 
semiregulars (2.1 x 10~® and 4.6 x 10“^ Mq yr“^ respec¬ 
tively), the separation is not as marked as in the color-period 
space. A third possibility is that the separation in color re¬ 
flects different circumstellar dust temperatures. Given that 
the IRAC bands are sensitive to dust in the temperature 
range 500 — 800 K, the observed dichotomy may indicate 
that semiregulars are deficient in the hot dust emission that 
is present in Mira variables. Thi s is consistent with the sug¬ 
gestion in iMarengo et all |2003) that semiregular variables 
could have discontinuous mass-loss, spending a significant 
part of their time in the quiescent phase with only brief 
phases of higher mass-loss. According to this hypothesis, 
semiregulars would have more excess during these brief in¬ 
tervals of higher mass-loss, but once this phase has con¬ 
cluded, the dust would expand and cool, leading to a small 
[3.6] —[8.0] excess (< 0.5 mag). Infrared color variations due 
to intermittent dust production have been found for many 
AGB stars (see e.g. IWhitelock et al.l l2006l . and references 
therein), for both Mira and semiregular LPVs. Those vari¬ 
ations have been attributed to the emission of dust clumps, 
similar to the ones observed with interferometric observa- 
tions of the very obscured Mira carbon star IRC-l-10216 
dTuthill et al.l[200ol ~). The different color we observe between 
Mira and SR variables may reflect the fact that dust episodes 
in semiregulars happen on longer time-scales, giving the dust 
time to cool between each mass-loss episode. Less excess for 
the semiregulars would mean that their IRAC colors are a 
better reflection of the star’s intrinsic color (and not thermal 
emission from the dust). This would lead to the more strict 
period-color relation for semiregular than for Mira variables. 
Indeed it is the higher excess in Mira variables that would 
prevent a similarly tight relation. 


6 MASS LOSS AND COMPARISON TO 
LMC/SMC AGB STARS 

The intense mass-loss processes that are often ac¬ 
tive during the AGB phase are not well under¬ 
stood. Co mmonly used e mpirical mass-l o ss-rate formulas 
(see e.g. Sal peterl 1974: Rei mersI 19751: Baud fc Habing 


1983l:lNieuwenhuiizen fc De Jagerj|l990l:IVassiliadis fc Wood 


1993f ) only provide an order-of-magnitude estimate of the 
mass-loss rate of AGB stars and cannot predict the actual 
mass-loss rate as a function of stellar parameters. Radiative 
transfer modeling can be used to infer the amount of cir¬ 
cumstellar dust from excess flux measurements in the ther¬ 
mal infrared. Coupled with assumptions about the AGB 
wind velocity (usually ~ 10 — 20 km s“^) and the gas-to- 
dust mass ratio (typically in the range of the ISM ratio of 
~ 200), these models can be used to provide estimates for 
AGB mass-loss rates from infrared photometric measure- 
me nts. This can be achieved either by fitting the SED (see 
e.g. lGullieuszik et al.ll201^ : lRiebel et al.l[2012f ) of each source 
to a model or adopting a color - mass-loss rel ation (see e.g. 
models using realistic stellar atmospheres from Groenewegen] 
I 2 OO 6 I : iGroenewegen et ^ l2009l . l2007l : iGruendl et al.ll2008ll . 
These models show how high-mass-loss rate AGB stars of 
all types follow monotonically increasing sequences similar 


to those shown in Figure [2] (higher mass-loss rate, stronger 
infrared excess). 

Figure fTT] shows the mass-loss rate plotted as function 
of [3.6] —[8.0] color. Reflecting the color separation found 
in Figure 1101 all semiregular variables are grouped in the 
0 — 0.5 mag color range. However, the overall mass-loss rates 
are comparable for the two groups. This is in accordance 
with the previously mentioned results found in Marengo et 
al. (2001; see discussion in Section [5l), and since we would 
expect semiregular variables to show less excess than Mira 
variables if they are characterised by discontinuous mass- 
loss. 

The overlap in the distribution of AGB stars of different 
chemical types echoes our finding that the majority of O- 
rich AGB stars in our sample have colors consistent with 
carbon stars (see Section 2)). Because our sample of O-rich 
AGB stars is biased towards higher mass-loss rates, the most 
significant contributors of O-rich AGB dust are also the most 
likely to be mis-classified with color cuts alone, suggesting 
that care must be taken in inferring the chemical yield of 
AGB dust from color-selected stellar population studies. 

Selecting for AGB stars with measured mass-loss rates 
may also explain the relatively high mass-loss rate that is 
characteristic of our sample. Our best-fit for the mass-loss 
rate as a function of [3.6] —[8.0] color for the O-rich AGB 
(excluding the supergiants) and C-rich AGB samples respec¬ 
tively, are: 

G — rich : 

log(dM/dt) = -19.20/[([3.6] - [8.0]) -k 3.65] - 1.45 , . 

O — rich : 

log(dM/dt) = -9.65/[([3.6[ - [8.0]) -k 1.35] - 2.29 

While the small number of sources in either sample does 
not allow us to establish meaningful uncertainties to the 
fits, we note that the two curves approach each other for 
[3.6]-[8.0] > 2 mag. This suggests that the mass-loss rate for 
red sources is similar, for a given value of the infrared excess, 
regardless of the dust chemistry. Our three supergiants with 
reliable mass-loss rate estimates follow a trend similar to 
the AGB stars, but as expected are shifted towards higher 
mass-loss rates. 

Similar fits are provided by different authors for AGB 
stars in the Magellanic Gloud s . Figure IIII shows a selection 
of them jMatsuu ra et al.ll2009l . I 2 OI 3 I : iGullieuszik et al.ll2012l : 

iRiebel et al.ll2012l ). based on the photometry obtained as 

part of the Spitzer/YRAC SAGE program. It is worth noting 
that the spread between the individual fits is large, as much 
as an order of magnitude in dM/dt. This spread may be, in 
part, a consequence of the specific selection criteria adopted 
by different authors, leading to different biases in the AGB 
samples used in the fit. Different assumptions for the gas- 
to-dust ratio and wind velocity, as well as different optical 
constants used to model the sources will also contribute to 
the spread. While our best-fitting relations predict larger 
mass-loss rates for our Galactic AGB sample, the large un¬ 
certainty in the fit parameters and the spread between the 
Magellanic Cloud curves prevent us from deriving meaning¬ 
ful conclusions about the dependence of mass-loss rates on 
metallicity. Such a dependence is expected for O-rich AGB 
stars, but not for C-rich sources, since carbon is synthesized 
locally in the Therma lly-Pulsing AGB stars, as argued by 
iMatsuura et al.l (l2009l L However, it is unclear whether the 
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Figure 11. A comparison of fits to the mass-loss rate as a func¬ 
tion of [3.6] —[8.0] color. The best-fitting relations to our sam¬ 
ple of Galactic C- and 0-rich AGB stars are shown with bold 
straight and dashed lines respectively. These are compared to fits 
found by other authors for larger samples of AGB stars in the 
LMC. Semiregular and irregular LPVs fall within the shaded re¬ 
gion while the Mira variables all have redder color. 

loc al synthesis of carbon increases the dust-to-gas ratio (see, 
e.g. lvan Loon et al.ll2008|) . A large Galactic sample with reli¬ 
able distances (that could be provided in the near future by 
GAIA) and radio-based mass-loss rate determinations will 

help resolve this issue. _ _ 

_ The selection criteria used by iMatsuura et al.l l|2009l. 

I2ni3h to differentiate evolved stars in the Magellanic Clouds 
result in a sample that is more directly comparable with 
our own Galactic sample. The sources in these works are 
identified by optical and infrared spectroscopy, avoiding the 
potential mis-classification of red 0-rich sources (OH/IR 
stars and RSG) in the extreme AGB class (where most 
of th e sources, at lea s t in the Magellanic Clouds, are C- 
rich). IMatsuura et al.l (l2013fl show how spectrally classified 
sources of different type tend to separate into three differ¬ 
ent regions in the [3.6] —[8.0] versus [8.0] color-magnitude 
diagram. These three regions are overlaid on our sample in 
Figure 1121 Note that for our Galactic sources, the bright 
supergiants also have the highest mass-loss rate plotted 
in Figure [TT] Our two 0-rich AGB stars with larger in- 
frared excess a r e in t he lowest luminosity region in the 
IMatsuura et al.l ll2013ll diagram, within the region where 
extreme C-rich AGB stars are expected. However, these 
two stars are characterised by a mass-loss rate as high as 
C-rich AGB stars wit h sim ilar excess. As poin ted out by 
IMatsuura et al.l ll2013ll and iRiebel et al.l (I 2 OI 2 I '). the over¬ 
all dust return of AGB stars to the fSM of the Magellanic 
Clouds is dominated by a handful of stars with very large 
mass-loss rates. In the Magellanic Clouds, most of these 
stars tend to be C-rich ‘extreme’ stars. The situation may 
be different in the Galaxy where a large number of very red 
OH/IR stars are found. These O-rich evolved stars, with 
mass-loss rates as high as the two reddest O-rich AGB stars 
in our sample, may contribute as much as the C-rich AGB 
stars to the overall dust budget in our Galaxy. A larger por- 



Figure 12. Color-Magnitude diagram of IRAC band [8.0]. Aver¬ 
age error is plotted in bot t om c entre. Color-magnitude selection 
cuts from IMatsuura et al.l ||2009| 'I overplotted. 

tion of O-rich AGB stars with high mass-loss rates in sources 
with near solar met allicity has indeed been observed in M33 
l|.Iavadi et al.l[ 201 ^ . 


7 CONCLUSIONS 

We present the results of the first study to characterise 
Galactic AGB stars in the IRAG bands. The sample con¬ 
sists of 48 AGB stars observed in two epochs - 22 O-rich, 
19 carbon-rich, and 7 S-type - representing a diverse AGB 
population with well-determined distance, chemistry, vari¬ 
ability type, and mass-loss rate. These are compared to 6 
supergiants representing mass-losing evolved stars outside 
the AGB. 

By fitting a PSF to unsaturated parts of the IRAC im¬ 
ages, we derive reliable photometry in all four IRAC bands. 
From this, we compute the mid-IR colors of Galactic AGB 
stars. AGB stars with O-rich chemistry (including S-type 
AGB stars) tend to have redder [3.6] —[8.0] colors than car¬ 
bon stars for a given [3.6] —[4.5] color, possibly driven by 
silicate emission in the [8.0] band. For colors including the 
5.8 [rm band, carbon stars separate into two distinct se¬ 
quences. Garbon stars with higher mass-loss rates tend to 
lie along the redder branch, while sources along the blue 
branch have lower mass-loss rates. We interpret this as evi¬ 
dence of a photospheric C 3 feature that is only visible in the 
low-mass-loss-rate sources that are relatively unobscured by 
dusty circumstellar envelopes. 

AGB stars are LPVs, so we have examined both the 
color and the magnitude as a function of period. We find a 
period-color re l ations hip consistent with the hypothesis of 
iMarengo et al.l (I 2 OOII I that semiregular variables lose mass 
discontinuously, leading to a lack of hot dust producing red¬ 
der colors. In contrast, the Mira variables extend along a 
more linear sequence toward redder colors as we might ex¬ 
pect based on the fact that their mass-loss rates are more 
sustained in time. 

The period-luminosity distribution of Galactic AGB 
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stars is similar to that found bv iRiebel et alJ ll2012ll for the 
LMC. Mira variables fall along the fundamental pulsation 
sequence C. Semiregular variables are mostly in sequences 
A and B, reflecting the presence of higher overtone modes 
and shorter periods. 

We derive a relationship between mass-loss rate and 
[3.6] —[8.0] color separately for 0-rich AGB stars and carbon 
stars. The overall shape of the fits is similar to those found 
by other authors for AGB stars in the LMG, albeit corre¬ 
sponding to a higher mass-loss rate for a given [3.6] —[8.0] 
color. Discrepancies between our fits and those found by 
other authors likely reflect different assumptions used in the 
radiative transfer models used to derive mass-loss rates and 
uncertainties in the gas-to-dust ratio. In addition, we find 
that neither color nor mass-loss rate is a good discriminator 
of chemical type, suggesting that care must be taken when 
inferring the chemical contribution of dust returned to the 
ISM by AGB stars. 
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Table 1. IRAC AGB OBSERVED TARGETS 


Target 

IRAS name 

Var. 

Period 

[days] 

Dist 

[kpc] 

Dist 

ref. 

M 

[Mo/yr] 

M 

ref. 

J Mag 

H Mag 

K Mag 

O-RICH STARS 

KU And 

00042+4248 

M 

750 

0.68 

14 

2.4x10-5 

11 

3.041+0.216 

1.829+0.156 

1.115+0.208 

RW And 

00445+3224 

M 

430 

0.86 

7 

1.9x10-’^ 

17 

3.052+0.238 

2.225+0.210 

1.765+0.228 

VY Gas 

00484+6238 

SRb 

100 

0.57 

16 



2.471+0.254 

1.484+0.186 

1.114+0.234 

SV Psc 

01438+1850 

SRb 

102 

0.38 

12 

3.0xl0-'^t 

12 

2.016+0.244 

1.012+0.190 

0.722+0.188 

RR Per 

02251+5102 

M 

389 

0.82 

15 



3.217+0.244 

2.325+0.196 

1.662+0.244 

RV Gam 

04265+5718 

SRb 

101 

0.35 

12 

2.5x10"'^ 

12 

1.667+0.254 

0.602+0.164 

0.412+0.170 

ET Vir 

14081-1604 

SRb 

80 

0.14 

16 



1.651+0.266 

0.814+0.290 

0.598+0.304 

RW Boo 

14390+3147 

SRb 

209 

0.29 

16 



1.560+0.206 

0.530+0.150 

0.301+0.170 

AX Sco 

16387-2700 

SRb 

138 

0.34 

16 



2.159+0.270 

1.224+0.282 

0.901+0.320 

TV Dra 

17081+6422 

SR 


0.51 

16 

2 .0x10-’^ 

11 

1.942+0.416 

1.058+0.174 

0.693+0.194 

V438 Oph 

17123+1107 

SRb 

169 

0.42 

4 



1.871+0.244 

0.935+0.292 

0.551+0.312 

TY Dra 

17361+5746 

Lb 


0.43 

12 

6 .0x10“’^ 

12 

2.410+0.296 

1.465+0.158 

1.083+0.194 

GZ Ser 

18347-0241 

Lb 


0.44 

12 

8 .0x10-’^ 

12 

2.557+0.302 

1.405+0.274 

0.930+0.310 

FI Lyr 

18401+2854 

SRb 

146 

0.88 

6 

3.8x10-® 

6 

2.040+0.282 

1.121+0.164 

0.715+0.212 

V1351 Gyg 

19409+5520 

Lb 


0.32 

16 



1.816+0.266 

0.787+0.170 

0.561+0.208 

Z Gyg 

20000+4954 

M 

263 

0.94 

5 

4.0x10-® 

19 

4.176+0.268 

3.258+0.228 

2.557+0.306 

V1300 Aql 

20077-0625 

M 

680 

0.66 

11 

1.3x10-5 

11 

6.906+0.032 

3.923+0.260 

2.059+0.262 

V584 Aql 

20079-0146 

Lb 


0.39 

12 

5.0x10-® 

12 

2.136+0.312 

1.147+0.204 

0.801+0.266 

RX Vul 

20507+2310 

M 

457 

0.71 

18 



2.724+0.290 

1.655+0.188 

1.083+0.228 

UX Gyg 

20529+3013 

M 

565 

1.85 

9 

3.2x10-5 

11 

2.910+0.310 

1.887+0.180 

1.400+0.194 

SS Peg 

22315+2418 

M 

424 

0.71 

18 



2.549+0.270 

1.513+0.182 

1.121+0.162 

V563 Gas 

23147+6009 

M 

534 

2.09 

11 



5.088+0.037 

3.882+0.036 

3.227+0.248 

C-RICH STARS 

VX And 

00172+4425 

SRa 

369 

0.56 

6 

1.4x10-^ 

6 

3.187+0.266 

1.891+0.202 

1.193+0.202 

HV Gas 

01080+5327 

M 

527 

0.97 

14 

3.2x10-5 

6 

5.585+0.039 

3.911+0.222 

2.466+0.290 

R For 

02270-2619 

M 

339 

0.97 

14 

5.0x10-5 

6 

4.230+0.274 

2.537+0.206 

1.349+0.274 

V623 Gas 

03075+5742 

Lb 


0.51 

6 

l.lxlO-’’ 

6 

2.869+0.210 

1.706+0.204 

1.208+0.212 

UV Aur 

05185+3227 

M 

394 

1.09 

6 

1.6x10-5 

6 

4.029+0.214 

3.018+0.192 

2.129+0.220 

SY Per 

04127+4030 

SR 

474 

1.43 

6 

1.5x10-5 

6 

4.596+0.316 

3.033+0.274 

2.132+0.340 

TU Tan 

05421+2424 

SRb 

190 

1.05 

6 

5.7x10-'^ 

6 

3.331+0.266 

2.093+0.210 

1.574+0.220 

BN Mon 

06192+0722 

SRb 

600 

1.28 

1 

4.2x10-5 

6 

4.517+0.312 

3.096+0.256 

2.322+0.308 

GR Gem 

06315+1606 

Lb 


0.92 

6 

8.4x10-^ 

6 

3.575+0.302 

2.030+0.250 

1.538+0.240 

V614 Mon 

06585-0310 

SRb 

60 

0.48 

6 

1.5x10-® 

6 

3.227+0.206 

2.307+0.206 

1.764+0.224 

GGGS 6296 

08305-3314 



4.21 

6 

l.OxlO-"' 

6 



11.1 

SZ Gar 

09582-5958 

SRb 

126 

0.37 

6 

4.6x10-’' 

6 

3.206+0.332 

2.112+0.280 

1.600+0.346 

V GrB 

15477+3943 

M 

358 

0.84 

6 

1.3x10-5 

6 

3.474+0.272 

2.209+0.218 

1.321+0.276 

SX Sco 

17441-3541 

SR 


0.83 

6 

3.8x10-^ 

6 

3.104+0.222 

1.884+0.210 

1.390+0.254 

EX Ser 

18040-0941 

Lb 


1.23 

6 

3.1x10-5 

6 

7.136+0.023 

4.720+0.076 

2.702+0.218 

DR Ser 

18448+0523 

Lb 


1.29 

6 

1.0x10-5 

6 

3.833+0.270 

2.488+0.218 

1.839+0.286 

S Set 

18476-0758 

SRb 

148 

0.58 

6 

5.6x10-5 

6 

2.303+0.314 

1.140+0.262 

0.627+0.288 

RS Gyg 

20115+3834 

SRa 

418 

0.65 

6 

2 .0x10-’’ 

6 

3.206+0.228 

1.985+0.190 

1.321+0.246 

T Ind 

21168-4514 

SRb 

320 

0.65 

6 

1.7x10-’ 

6 

1.981+0.248 

0.970+0.262 

0.564+0.316 

S-RIGH STARS 

R Gem 

07043+2246 

M 

370 

0.71 

13 

2.1x10-5 

11 

2.530+0.244 

1.642+0.218 

1.459+0.226 

NQ Pup 

07507-1129 

Lb 


0.81 

16 



3.542+0.292 

2.587+0.212 

2.306+0.274 

S UMa 

12417+6121 

M 

226 

1.09 

18 



4.458+0.214 

3.431+0.196 

3.017+0.246 

RCyg 

19354+5005 

M 

426 

0.44 

14 

1 .6x10-’’ 

11 

2.251+0.314 

1.379+0.208 

0.861+0.246 

AA Gyg 

20026+3640 

SRb 

213 

0.48 

13 



2.066+0.288 

1.050+0.166 

0.625+0.224 

X Aqr 

22159-2109 

M 

312 

1.32 

18 



4.710+0.190 

3.823+0.196 

3.338+0.248 

HR Peg 

22521+1640 

SRb 

50 

0.42 

16 



2.306+0.270 

1.239+0.194 

1.041+0.206 

SUPERGIANTS 

XX Per 

01597+5459 

SRc 

415 

2.29 

3 



3.442+0.262 

2.480+0.232 

1.972+0.242 

W Per 

02469+5646 

SRc 

485 

2.29 

3 



3.095+0.202 

1.999+0.172 

1.568+0.252 

NO Aur 

05374+3153 

Lc 


0.60 

16 



2.122+0.322 

1.128+0.202 

0.971+0.196 

U Lac 

22456+5453 

SRc 


2.75 

10 



3.815+0.224 

2.672+0.166 

1.932+0.206 

NML Gyg* 

R20445+3955 

SRc 

940t 

1.61 

20 

6.4x10-5 

21,22 

4.877+0.037 

2.389+0.200 

0.791+0.204 

VX Sgr* 

18050-2213 

SRc 

732 

1.57 

2 

6.1x10-5 

4 

1.744+0.260 

0.550+0.304 

-0.122+0.362 


* from Table 4.4 in^chust^ (^007), f multiple velocity components detected, | ^onnie^^al^^97. 

References: (1) ^ergea^^^Chevallier ^005; (2) ^hei^^al^^OOT; (3) ^urri^e^al^^O^; (4) 

(f\\ niior^/4clin; ol onnfi. /'7^ niior./4cUn5 Sr Riioor. OHOft- (Uii rr, i-^V. r-^-.ro IQTS- /'Q^ T^iir-o-..o rv, o cl '*5bnCVTl OV ol OntTC-' ' 
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Table 2. IRAC PHOTOMETRY 


Target 

IRAS name 

epoch 

MJD“ 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

O-RICH STARS 

KU And 

00042+4248 

1 

53961.9209 

1.140+0.031 

0.291+0.021 

-0.400+0.015 

-1.584+0.076 

KU And 

00042+4248 

2 

54328.9385 

1.747+0.027 

0.879+0.024 

0.108+0.018 

-1.175+0.074 

RW And 

00445+3224 

1 

53960.9362 

1.109+0.030 

1.042+0.028 

0.742+0.022 

0.216+0.066 

RW And 

00445+3224 

2 

54324.8287 

0.904+0.025 

0.831+0.023 

0.546+0.027 

0.096+0.059 

VY Gas 

00484+6238 

1 

54005.9215 

0.879+0.024 

0.981+0.027 

0.786+0.022 

0.564+0.027 

VY Gas 

00484+6238 

2 

54149.6210 

0.843+0.024 

0.929+0.026 

0.731+0.032 

0.485+0.025 

SV Psc 

01438+1850 

1 

53961.9248 

0.395+0.023 

0.468+0.025 

0.256+0.021 

-0.027+0.042 

SV Psc 

01438+1850 

2 

54149.6256 

0.305+0.015 

0.411+0.016 

0 .202+0.020 

-0.063+0.031 

RR Per 

02251+5102 

1 

54005.9269 

0.427+0.024 

0.320+0.015 

0.044+0.017 

-0.289+0.033 

RR Per 

02251+5102 

2 

54149.9993 

0.764+0.022 

0.775+0.022 

0.502+0.017 

0 .202+0.020 

RV Gam 

04265+5718 

1 

54005.9347 

0.327+0.015 

0.451+0.016 

0.277+0.021 

0.073+0.023 

RV Gam 

04265+5718 

2 

54396.8337 

0.284+0.014 

0.435+0.016 

0.216+0.020 

0.073+0.023 

ET Vir 

14081-1604 

1 

53957.3504 

0.395+0.016 

0.573+0.018 

0.443+0.024 

0.372+0.031 

ET Vir 

14081-1604 

2 

54148.6245 

0.387+0.023 

0.537+0.027 

0.427+0.024 

0.372+0.031 

RW Boo 

14390+3147 

1 

53922.4493 

0.127+0.018 

0.387+0.016 

0.145+0.019 

-0.053+0.032 

RW Boo 

14390+3147 

2 

54149.9923 

0.133+0.018 

0.387+0.016 

0.170+0.019 

-0.048+0.033 

AX Sco 

16387-2700 

1 

53997.8064 

0.564+0.018 

0.649+0.020 

0.435+0.016 

0.256+0.021 

AX Sco 

16387-2700 

2 

54357.9155 

0.582+0.019 

0.700+0.015 

0.485+0.017 

0.291+0.028 

TV Dra 

17081+6422 

1 

54005.9137 

0.177+0.021 

0.263+0.021 

0.061+0.017 

-0.123+0.019 

TV Dra 

17081+6422 

2 

54286.3724 

0.170+0.019 

0.242+0.020 

0.022+0.017 

-0.161+0.019 

V438 Oph 

17123+1107 

1 

54002.9634 

0.209+0.020 

0.294+0.014 

0.096+0.018 

-0.094+0.020 

V438 Oph 

17123+1107 

2 

54357.9197 

0.164+0.019 

0.263+0.021 

0.073+0.017 

-0.180+0.020 

TY Dra 

17361+5746 

1 

54005.9110 

0.679+0.020 

0.891+0.020 

0.649+0.030 

0.284+0.042 

TY Dra 

17361+5746 

2 

54286.3697 

0.689+0.020 

0.879+0.024 

0.659+0.030 

0.298+0.057 

GZ Ser 

18347-0241 

1 

54005.8978 

0.519+0.026 

0.710+0.021 

0.435+0.024 

0.277+0.028 

GZ Ser 

18347-0241 

2 

54229.6547 

0.591+0.037 

0.742+0.032 

0.476+0.025 

0.291+0.028 

FI Lyr 

18401+2854 

1 

54005.9018 

0.411+0.016 

0.601+0.028 

0.357+0.023 

0.108+0.030 

FI Lyr 

18401+2854 

2 

54357.9236 

0.419+0.024 

0.601+0.019 

0.372+0.023 

0.114+0.048 

V1351 Gyg 

19409+5520 

1 

53956.9174 

0.196+0.020 

0.395+0.023 

0.229+0.020 

0.177+0.019 

V1351 Gyg 

19409+5520 

2 

54286.3181 

0.209+0.020 

0.403+0.016 

0 .222+0.020 

0.164+0.025 

Z Cyg 

20000+4954 

1 

54005.9076 

2.174+0.056 

1.892+0.043 

1.620+0.047 

0.929+0.077 

Z Gyg 

20000+4954 

2 

54286.3642 

1.862+0.060 

1.669+0.051 

1.422+0.060 

0.753+0.065 

V1300 Aql 

20077-0625 

1 

54430.9969 

0.145+0.025 

-0.483+0.028 

-1.187+0.022 

-2.280+0.040 

V1300 Aql 

20077-0625 

2 

54637.7548 

0.649+0.020 

0.050+0.017 

-0.647+0.024 

-1.627+0.036 

V584 Aql 

20079-0146 

1 

54064.4277 

0.427+0.016 

0.502+0.017 

0.305+0.022 

0.091+0.024 

V584 Aql 

20079-0146 

2 

54430.9864 

0.427+0.016 

0.502+0.017 

0.294+0.021 

0.033+0.034 

RX Vul 

20507+2310 

1 

54286.3354 

0.419+0.024 

0.468+0.017 

0.170+0.019 

-0.104+0.020 

RX Vul 

20507+2310 

2 

54430.9901 

0.493+0.026 

0.610+0.029 

0.342+0.022 

0.044+0.034 

UX Gyg 

20529+3013 

1 

54286.3327 

1.009+0.027 

0.731+0.021 

0.411+0.024 

-0.403+0.075 

UX Gyg 

20529+3013 

2 

54430.9927 

1.629+0.034 

1.307+0.036 

0.968+0.026 

0.183+0.051 

SS Peg 

22315+2418 

1 

54286.3777 

0.820+0.023 

0.649+0.020 

0.380+0.023 

0.005+0.022 

SS Peg 

22315+2418 

2 

54100.2776 

0.669+0.030 

0.669+0.020 

0.419+0.016 

0.164+0.038 

V563 Gas 

23147+6009 

1 

53956.1633 

2.585+0.035 

2.097+0.037 

1.684+0.041 

1.036+0.056 

V563 Gas 

23147+6009 

2 

54149.4121 

2.935+0.032 

2.639+0.037 

2.190+0.041 

1.620+0.097 

G-RIGH STARS 

VX And 

00172+4425 

1 

53961.0086 

0.158+0.019 

0.357+0.023 

0.485+0.025 

-0.216+0.027 

VX And 

00172+4425 

2 

54150.6613 

0.183+0.019 

0.395+0.023 

0.528+0.018 

-0.157+0.019 

HV Gas 

01080+5327 

1 

53961.4766 

0.891+0.025 

0.342+0.015 

-0.108+0.029 

-0.589+0.032 

HV Gas 

01080+5327 

2 

54150.6238 

0.357+0.015 

-0.152+0.019 

-0.573+0.016 

-0.955+0.023 

R For 

02270-2619 

1 

53961.0145 

-0.476+0.018 

-0.900+0.019 

-1.227+0.018 

-1.621+0.037 

R For 

02270-2619 

2 

54357.9070 

-0.354+0.016 

-0.811+0.016 

-1.160+0.019 

-1.571+0.051 

V623 Gas 

03075+5742 

1 

54005.9318 

0.485+0.017 

0.546+0.018 

0.305+0.022 

0.102+0.018 

V623 Gas 

03075+5742 

2 

54357.9005 

0.493+0.017 

0.564+0.018 

0.357+0.023 

0.120+0.018 

SY Per 

04127+4030 

1 

54396.8364 

0.820+0.046 

0.786+0.034 

0.649+0.040 

0.605+0.029 

SY Per 

04127+4030 

2 

54531.1417 

0.942+0.039 

0.879+0.037 

0.775+0.033 

0.387+0.047 

UV Aur 

05185+3227 

1 

54005.9383 

1.422+0.029 

1.017+0.022 

0.620+0.019 

0.164+0.025 

UV Aur 

05185+3227 

2 

54191.1691 

1.171+0.026 

0.855+0.024 

0.451+0.025 

0.170+0.025 

TU Tau 

05421+2424 

1 

54396.8428 

0.981+0.027 

1.094+0.029 

0.929+0.038 

0.411+0.024 

TU Tau 

05421+2424 

2 

54188.7098 

0.955+0.021 

1.051+0.029 

0.916+0.025 

0.395+0.023 

BN Mon 

06192+0722 

1 

54038.7664 

1.463+0.029 

1.550+0.036 

1.422+0.040 

1.036+0.042 

BN Mon 

06192+0722 

2 

54396.8610 

1.454+0.021 

1.550+0.032 

1.463+0.042 

1.022+0.055 
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Table 2 — continued 


Target 

IRAS name 

epoch 

MJD“ 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

C-RICH STARS 

CR Gem 

06315+1606 

1 

54061.8204 

0.995+0.027 

0.916+0.025 

0.601+0.028 

0.485+0.034 

CR Gem 

06315+1606 

2 

54396.8582 

0.855+0.024 

0.867+0.024 

0.555+0.027 

0.419+0.032 

V614 Mon 

06585-0310 

1 

54069.0020 

1.272+0.025 

1.422+0.028 

1.254+0.026 

0.955+0.039 

V614 Mon 

06585-0310 

2 

54190.6094 

1.272+0.025 

1.422+0.032 

1.265+0.028 

1.009+0.041 

CGCS 6296 

08305-3314 

1 

54466.0597 

4.120+0.048 

2.447+0.052 

1.124+0.031 

-0.133+0.029 

CGCS 6296 

08305-3314 

2 

54600.7407 

4.674+0.121 

2.935+0.049 

1.550+0.045 

0.183+0.064 

SZ Gar 

09582-5958 

1 

53888.8675 

1.065+0.029 

1.171+0.032 

1.065+0.029 

0.742+0.032 

SZ Gar 

09582-5958 

2 

54148.6298 

0.976+0.027 

1.080+0.029 

0.916+0.038 

0.659+0.030 

V CrB 

15477+3943 

1 

53959.8338 

0.874+0.019 

0.387+0.016 

-0.048+0.021 

-0.385+0.015 

V CrB 

15477+3943 

2 

54152.9619 

0.403+0.016 

-0.063+0.015 

-0.447+0.014 

-0.811+0.021 

SX Sco 

17441-3541 

1 

53997.2559 

0.891+0.049 

1.022+0.042 

0.942+0.039 

0.476+0.034 

SX Sco 

17441-3541 

2 

54396.9331 

0.867+0.048 

1.065+0.043 

0.995+0.041 

0.476+0.034 

FX Ser 

18040-0941 

1 

54002.9673 

-0.118+0.019 

-0.787+0.016 

-1.298+0.023 

-1.811+0.041 

FX Ser 

18040-0941 

2 

54396.9368 

0.263+0.021 

-0.433+0.022 

-0.969+0.022 

-1.505+0.041 

DR Ser 

18448+0523 

1 

54039.4356 

1.344+0.037 

1.505+0.043 

1.402+0.040 

1.051+0.043 

DR Ser 

18448+0523 

2 

54229.1684 

1.442+0.041 

1.572+0.046 

1.505+0.065 

1.156+0.063 

S Set 

18476-0758 

1 

54005.8936 

-0.058+0.021 

0.096+0.024 

0.044+0.028 

-0.338+0.020 

S Set 

18476-0758 

2 

54396.9396 

-0.099+0.020 

0.073+0.023 

0.028+0.028 

-0.385+0.038 

RS Cyg 

20115+3834 

1 

54036.9492 

0.963+0.026 

1.080+0.029 

0.867+0.024 

0.764+0.055 

RS Cyg 

20115+3834 

2 

54286.3299 

1.254+0.034 

1.065+0.043 

0.742+0.022 

0.867+0.048 

T Ind 

21168-4514 

1 

54069.0284 

0.242+0.020 

0.380+0.015 

0.177+0.019 

0.005+0.022 

T Ind 

21168-4514 

2 

54430.9820 

0.263+0.021 

0.387+0.023 

0.216+0.020 

0.028+0.017 

S-RICH STARS 

R Gem 

07043+2246 

1 

54068.8159 

1.210+0.027 

1.214+0.027 

0.891+0.037 

0.476+0.051 

R Gem 

07043+2246 

2 

54228.7301 

1.140+0.031 

1.080+0.029 

0.786+0.022 

0.468+0.033 

NQ Pup 

07507-1129 

1 

54069.3910 

1.990+0.034 

2.190+0.041 

2.024+0.065 

1.990+0.034 

NQ Pup 

07507-1129 

2 

54228.8126 

1.990+0.034 

2.215+0.042 

2.024+0.070 

1.990+0.068 

S UMa 

12417+6121 

1 

54069.7437 

2.620+0.037 

2.834+0.030 

2.627+0.024 

2.468+0.053 

S UMa 

12417+6121 

2 

53889.3814 

3.002+0.035 

3.189+0.041 

3.002+0.052 

2.798+0.072 

RCyg 

19354+5005 

1 

54005.9052 

0.067+0.017 

0.022+0.017 

-0.377+0.015 

-0.856+0.049 

RCyg 

19354+5005 

2 

54286.3666 

1.124+0.031 

1.009+0.041 

0.476+0.034 

-0.079+0.051 

AA Cyg 

20026+3640 

1 

54065.7245 

0.202+0.013 

0.528+0.018 

0.236+0.020 

0.139+0.019 

AA Cyg 

20026+3640 

2 

54286.3611 

0.216+0.020 

0.573+0.018 

0.256+0.021 

0.177+0.032 

X Aqr 

22159-2109 

1 

54286.3244 

2.190+0.065 

2.135+0.062 

1.911+0.044 

1.721+0.053 

X Aqr 

22159-2109 

2 

54062.9024 

2.729+0.047 

2.677+0.045 

2.447+0.052 

2.215+0.042 

HR Peg 

22521+1640 

1 

53928.4751 

0.649+0.020 

0.891+0.025 

0.710+0.031 

0.620+0.038 

HR Peg 

22521+1640 

2 

54095.8618 

0.639+0.020 

0.843+0.024 

0.669+0.030 

0.582+0.028 

SUPERGIANTS 

XX Per 

01597+5459 

1 

54005.9244 

1.171+0.032 

1.088+0.030 

0.721+0.032 

0.085+0.059 

XX Per 

01597+5459 

2 

54150.6211 

1.220+0.033 

1.156+0.032 

0.786+0.034 

0.145+0.062 

W Per 

02469+5646 

1 

54005.9295 

1.237+0.036 

1.171+0.032 

0.808+0.034 

0.177+0.064 

W Per 

02469+5646 

2 

54152.9557 

1.244+0.034 

1.178+0.032 

0.820+0.035 

0.114+0.060 

NO Aur 

05374+3153 

1 

54396.8402 

0.610+0.019 

0.843+0.024 

0.582+0.019 

0.320+0.029 

NO Aur 

05374+3153 

2 

54188.5900 

0.620+0.019 

0.843+0.024 

0.591+0.019 

0.312+0.022 

U Lac 

22456+5453 

1 

54005.9181 

0.786+0.022 

0.573+0.018 

0.196+0.020 

-0.462+0.043 

U Lac 

22456+5453 

2 

54328.9353 

0.764+0.022 

0.564+0.018 

0.183+0.019 

-0.440+0.024 

NML Cyg* 

R20445+3955 

1 

53213.2738 

-1.53+0.15 

-2.42+0.15 


-4.53+0.22 

VX Sgr* 

18050-2213 

1 

53634.9152 

-1.26+0.15 

-1.83+0.15 


-3.36+0.15 


“ MDJ = JD-2,400,000.5 

* only one epoch; photometry from^chust^ ^007, see Table 4.4) 
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Table 3. Mira/SR P-L sequence 


Name 

Type 

Sequence 

Notes 

O-RICH STARS 

KU And 

M 

D 

low [3.6], large IR excess — extinction from the circumstellar envelope? 

RW And 

M 

C 

along the bottom of the sequence 

VY Cas 

SRb 

B-C’ 


SV Psc 

SRb 

C’ 


RR Per 

M 

C 


RV Cam 

SRb 

C’ 


ET Vir 

SRb 

C 

smallest 3.6 p.m magnitude 

RW Boo 

SRb 

C 

at the bottom of sequence C, small IR excess 

AX Sco 

SRb 

C 

near the top of sequence C 

V438 Oph 

SRb 

C’-C 


FI Lyr 

SRb 

B 

at the top of sequence B 

Z Cyg 

M 

C 

large IR excess 

V1300 Aql 

M 

C-D 

large IR excess, low 3.6 |i.m, extinction from the circumstellar envelope? 

RX Vul 

M 

C 

large IR excess, lies at the bottom of C 

UX Cyg 

M 

C 

one of largest IR excesses 

SS Peg 

M 

c+ 

just below sequence C, modest 3.6 \im and IR excess compared to other Miras of similar period 

V563 Cas 

M 

C-D 

just below C, among largest IR excesses 

C-RICH STARS 

VX And 

SRa 

C 


HV Cas 

M 

C 


R For 

M 

C’ 

highest Mira 3.6 pm magnitude 

UV Aur 

M 

C 


SY Per 

SR 

C 


TU Tau 

SRb 

C’ 

largest IR excess of the SRs 

BN Mon 

SRb 

C-D 

falls just below sequence C 

V614 Mon 

SRb 

B 


CGCS 6296 



no period available 

SZ Car 

SRb 

C 


V CrB 

M 

C 


SX Sco 

SR 


no period available 

S Set 

SRb 

B 


RS Cyg 

SRa 

C-D 

one of two SRs below C 

T Ind 

SRb 

C 


S-RICH STARS 

R Gem 

M 

C-D 

just below C, modest 3.6 pm and IR excess compared to other Miras of similar period 

S UMa 

M 

C 

at the top of sequence C, smallest IR excess of the Miras 

RCyg 

M 

C-D 

slightly closer to D than C 

AA Cyg 

SRb 

C 

at the top of sequence C 

X Aqr 

M 

C 

second smallest Mira IR excess 

HR Peg 

SRb 

A 

only source in A 

SUPERGIANTS 

XX Per 

SRc 

C’ 

SR with large IR excess 

W Per 

SRc 

C’ 

right along the lower edge of sequence C’; SR with large IR excess 

NML Cyg 

SRc 

C’ 

brightest supergiant 

VX Sgr 

SRc 

B-C’ 

along the boundary between sequences B and C’; second brightest supergiant 
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